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Abstract. In this study the effect of dust aerosol on up- 
per tropospheric cirrus clouds through heterogeneous ice 
nucleation is investigated in the Community Atmospheric 
Model version 5 (CAM5) with two ice nucleation param- 
eterizations. Both parameterizations consider homogeneous 
and heterogeneous nucleation and the competition between 
the two mechanisms in cirrus clouds, but differ significantly 
in the number concentration of heterogeneous ice nuclei 
(IN) from dust. Heterogeneous nucleation on dust aerosol re- 
duces the occurrence frequency of homogeneous nucleation 
and thus the ice crystal number concentration in the North- 
ern Hemisphere (NH) cirrus clouds compared to simulations 
with pure homogeneous nucleation. Global and annual mean 
shortwave and longwave cloud forcing are reduced by up 
to 2.0 zb 0.1 Wm -2 (la uncertainty) and 2.4 zb 0.1 Wm -2 , 
respectively due to the presence of dust IN, with the net 
cloud forcing change of —0.40 ± 0.20 Wm -2 . Comparison 
of model simulations with in situ aircraft data obtained in 
NH mid-latitudes suggests that homogeneous ice nucleation 
may play an important role in the ice nucleation at these re- 
gions with temperatures of 205-230 K. However, simulations 
overestimate observed ice crystal number concentrations in 
the tropical tropopause regions with temperatures of 190- 
205 K, and overestimate the frequency of occurrence of high 
ice crystal number concentration (> 200 L _1 ) and underes- 


timate the frequency of low ice crystal number concentra- 
tion (< 30 L -1 ) at NH mid-latitudes. These results highlight 
the importance of quantifying the number concentrations and 
properties of heterogeneous IN (including dust aerosol) in the 
upper troposphere from the global perspective. 


1 Introduction 

Cirrus clouds composed of ice crystals cover roughly 30 % 

of the globe (Wang et al., 1996; Wylie and Menzel, 1999), 

and have been shown to play an important role in modify- 
ing the global radiative balance through scattering shortwave 
radiation and absorbing and emitting longwave (LW) terres- 
trial radiation (Corti et al., 2005; Kay and Gettelman, 2009; 
Liou, 1986; Lohmann and Roeckner, 1995; Ramanathan and 

Collins, 1991). Cirrus clouds also modulate water vapor in 
the upper troposphere and lower stratosphere (Gettelman et 

al., 2002; Korolev and Isaac, 2006; Kramer et al., 2009). 

Although cirrus clouds are an important player in the 

global climate system, there are still large uncertainties in 
the understanding of cirrus cloud processes (DeMott et al., 

2003; Lohmann et al., 2008; Barahona and Nenes, 2011). 
Compared to droplet activation in warm clouds, ice nu- 
cleation in cold clouds is more complicated. There are 
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multiple ice nucleation mechanisms identified depending on 
aerosol surface properties and environmental conditions. Cir- 
rus clouds can be formed by the homogeneous freezing of 
cloud droplets and aerosol solution particles at low tempera- 
tures (< — 37 °C) and high relative humidities over ice (RHi, 
140-160%) (e.g., Koop et al., 2000). Cirrus clouds can also 
be formed by the heterogeneous nucleation of ice on aerosol 
particles (Pruppacher and Klett, 1997), called ice nuclei (IN). 

Various insoluble or partly insoluble aerosol particles can 
act as IN in deposition, contact, immersion, and condensation 
“freezing” modes, such as mineral dust, metallic particles, 
and primary biological particles (Cziczo et al., 2004; De- 
Mott et al., 2000; Szyrmer and Zawadzki, 1997). Although 
very scarce, heterogeneous IN could significantly impact the 
ice crystal number concentration compared to pure homo- 
geneous freezing by initiating ice nucleation earlier, deplet- 
ing available water vapor and ultimately hindering the occur- 
rence of homogeneous freezing (Chen et al., 2000; DeMott 
et al., 1994; Karcher and Lohmann, 2003; Liu and Penner, 
2005; Barahona and Nenes, 2009a, b). The relative impor- 
tance of heterogeneous IN depends on the concentration of 
IN as well as the cooling rate (DeMott et al., 1994; Karcher 
and Lohmann, 2003). Till now the importance of heteroge- 
neous IN on global upper tropospheric ice number concen- 
trations remains unclear, partially due to the scarcity of IN 
measurements and instrument artifacts of in situ ice crys- 
tal number measurements (Heymsfield, 2007; McFarquhar 
et al., 2007). Recent measurements indicate that heteroge- 
neous IN might play a more important role than was previ- 
ously thought in the Tropical Tropopause Layer (TTL) with 
extreme low temperatures (less than — 70°C) after correct- 
ing the measurement artifacts (Jensen et al., 2010) and after 
identifying new classes of particles and phase states that can 
act as IN (Abbatt et al., 2006; Murray et al., 2010). 

There have been significant advancements in the devel- 
opment of ice nucleation parameterizations for large-scale 
models. Homogeneous freezing of aqueous solutions can be 
relatively well represented by theory (Koop et al., 2000). For 
heterogeneous nucleation, empirical fits of IN observations 
are often used (Cooper, 1986; Cotton et al., 1986; DeMott et 
al., 1998; Meyers et al., 1992). Heterogeneous ice nucleation 
rates can also be calculated using the classical nucleation 
theory (CNT) (Fletcher, 1962; Pruppacher and Klett, 1997; 
Khvorostyanov and Curry, 2000, 2004). CNT is a statistical 
method that parameterizes the rate of heterogeneous ice nu- 
cleation as a function of ambient conditions such as tempera- 
ture and ice supers aturation, and properties of the IN such as 
size, contact angle of ice germ on the substrate and activation 
energy. Key parameters in CNT such as contact angle and 
activation energy can also be constrained using laboratory 
data for different types of aerosols (Chen et al., 2008), and 
possibly ambient observations (Curry and Khvorostyanov, 
2012). Phillips et al. (2008) derived empirical formulations 
for parameterizing atmospheric IN as a function of surface 
area densities of aerosol species (mineral dust, black car- 
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bon and hydrophobic organics) in addition to temperature 
and ice supersaturation, based on observations from several 
field campaigns obtained from the Colorado State Univer- 
sity (CSU) Count Flow Diffusion Chamber (CFDC). Using 
a larger dataset obtained from the same CFDC, DeMott et 
al. (2010) derived an empirical formulation of the number 
concentration of IN in mixed-phase clouds as a function of 
temperature and number concentration of aerosol particles 
with sizes larger than 0.5 pm. Ice crystal number concentra- 
tions predicted from the CNT are much higher than those 
from the empirical method by Phillips et al. (2008), when 
both methods were tested in the same cloud parcel model 
framework for the same assumed aerosol distribution (Eid- 
hammer et al., 2009). 

Determining the ice crystal concentration in cirrus clouds 
require (in addition to a nucleation theory) knowledge of the 
supersaturation that develops in clouds; this is typically ac- 
complished by considering a dynamical parcel framework 
(e.g., Pruppacher and Klett, 1997) and determining ice crys- 
tal number from the maximum supersaturation that develops. 
Karcher and Lohmann (2002a, b) first derived a formula- 
tion for homogeneous freezing of sulfate solution droplets 
based on this approach. The parameterization was extended 
in Karcher et al. (2006) to consider the competition of het- 
erogeneous ice nucleation on mineral dust and black carbon 
with homogeneous nucleation. Liu and Penner (2005) pre- 
sented a parameterization for ice number concentration for 
large-scale models based on parcel model numerical sim- 
ulations, and considered the competition between homoge- 
neous and heterogeneous nucleation at temperature less than 
— 35 °C. Heterogeneous nucleation was based on the CNT. 
The ice crystal concentration is related to updraft veloc- 
ity, temperature, and aerosol concentration. Barahona and 
Nenes (2009a, b) provided a framework for ice nucleation 
parameterizations by solving analytically the parcel model 
equations. These parameterizations provide expressions for 
ice crystal number concentration, explicitly considering the 
effects of aerosol size and number, updraft velocity, and de- 
position coefficient, and the competition between homoge- 
neous nucleation and heterogeneous nucleation. The IN nu- 
cleation “spectrum” (i.e., the number of IN as a function of 
parcel supersaturation) can have any forms and be derived 
from CNT or empirical data. 

Using the Liu and Penner (2005) parameterization imple- 
mented in the Community Atmospheric Model version 3 (Liu 
et al., 2007), Liu et al. (2009) investigated the impacts of an- 
thropogenic aerosol (black carbon and sulfate) on ice clouds 
and found that the effect of anthropogenic black carbon on 
cirrus clouds through ice nucleation depends critically on its 
ice nucleation efficiency assumed in the model and thus the 
relative importance of homogeneous and heterogeneous ice 
nucleation. Using the European Centre/Hamburg (ECHAM) 
general circulation model with the Karcher et al. (2006) ice 
nucleation parameterization, Hendricks et al. (2011) stud- 
ied the effects of heterogeneous IN (mineral dust and black 
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carbon) on cirrus clouds, found that the annual zonal mean 
cloud ice number concentration is reduced by up to 20 % in 
the tropics and 1-10 % in the midlatitudes due to the hetero- 
geneous IN. Significant reduction (up to 10%) in the mean 
ice water content are modeled, which leads to the reduction 
of annual and zonal mean water vapor mixing ratio of up to 
5 % at cirrus levels. 

In this study, the ice nucleation parameterization of Bara- 
hona and Nenes (2009b, hereafter as BN09) was imple- 
mented in the Community Atmosphere Model version 5 
(CAM5), which is the atmospheric component of the Com- 
munity Earth System Model version 1 (CESM1). Together 
with the existing Liu and Penner (2005, hereafter as LP05) 
ice nucleation parameterization in CAM5, we investigate the 
sensitivity of cirrus clouds to different formulations of het- 
erogeneous ice nucleation (through its competition with ho- 
mogeneous freezing). Modeled ice crystal number concen- 
tration and RHi are compared with available observational 
data to understand the ice nucleation mechanisms in the up- 
per troposphere. Sensitivity tests for pure homogeneous and 
pure heterogeneous nucleation were performed to examine 
the importance of heterogeneous IN from dust aerosol. The 
paper is organized as follows. Model description is presented 
in Sect. 2. Section 3 compares modeled ice crystal number 
concentrations from different parameterizations. Section 4 
compares model results with observational data. Section 5 
presents atmospheric impacts of dust IN and conclusions are 
in Sect. 6. 


2 Model description 

2.1 CAM5 

This study uses CAM5, the latest released version of CAM 
(Neale et al., 2010). Compared to CAM version 3 (Collins 
et al., 2006a, b), several important physics components have 
been updated. CAM5 includes a two-moment stratiform 
cloud microphysics scheme (Morrison and Gettelman, 2008, 
hereafter as MG08; Gettelman et al., 2008). The model ra- 
diation scheme has been updated to the Rapid Radiative 
Transfer Model for GCM described by Iacono et al. (2008). 
The deep convective scheme has been modified following 
Neale et al. (2008). The shallow convective scheme follows 
that of Park and Bretherton (2009). The moist boundary 
layer scheme is based on Bretherton and Park (2009). The 
cloud condensation parameterization uses a new liquid cloud 
macrophysics closure described by Park et al. (2012). The 
aerosol treatment in the model uses the modal aerosol mod- 
ule (MAM) (Liu et al., 2012). MAM treats the internal mix- 
ing within a mode and external mixing between modes of 
major aerosol species (sulfate, black carbon, primary organic 
carbon, secondary organic carbon, mineral dust and sea salt), 
and predicts mass mixing ratio of different aerosol species in 
each mode and number mixing ratio of each mode. Specifi- 
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cally, mineral dust important for this study is emitted in the 
accumulation mode with the size (diameter) range of 0.1- 
1 pm and coarse mode with the size (diameter) range of 1- 
10 pm following the emission scheme of Zender et al. (2003). 
Once in the atmosphere, mineral dust is transported horizon- 
tally and vertically, and is removed by the dry and wet depo- 
sition. On the annual average, ~ 62 % of the global emitted 
dust mass is removed by the dry deposition (mainly by the 
gravitational settling) predicted by MAM, while wet deposi- 
tion removes ~ 85 % of dust mass in the accumulation mode 
(Liu et al., 2012). 

2.2 Ice cloud scheme in CAM5 

Gettelman et al. (2010) modified cloud microphysics and 
macrophysics closure related to ice phase for CAM5 as com- 
pared to the CAM version 3.5 used in MG08. This includes: 
(1) condensation closure was modified to perform saturation 
adjustment only for liquid water. This allows the supersatu- 
ration over ice in ice clouds. Water vapor deposition onto ice 
crystals in mixed-phase and ice clouds is treated with pro- 
cess rate calculations. Sublimation occurs when water vapor 
is sub-saturated over ice; (2) the Bergeron-Lindeisen pro- 
cess is treated for mixed-phase clouds to deplete cloud liq- 
uid given the assumption that water vapor is saturated with 
respect to water in mixed-phase clouds; (3) the diagnostic 
fractional cloudiness scheme for ice clouds is still based on 
Slingo (1987), but calculated using the total water, not just 
water vapor; and (4) ice nucleation in MG08 is changed as 
discussed in Sect. 2.3. 

2.3 Ice nucleation parameterizations in the standard 
CAM5 

Along with vapor deposition, ice nucleation provides a 
source for cloud ice in MG08. Sink terms include aggrega- 
tion, autoconversion, accretion, melting and sedimentation. 
Aggregation reduces the ice particle number only and con- 
serves the ice water content, while the other sinks mentioned 
above affect both quantities. The empirical temperature- 
based IN parameterization from Cooper (1986) in MG08 was 
replaced in CAM5 with an explicit ice nucleation scheme 
with links to the aerosol concentration and type (Gettelman 
et al., 2010). Ice nucleation is based on Liu et al. (2007), 
which includes homogeneous nucleation on sulfate, hetero- 
geneous immersion freezing on mineral dust and competition 
between the two mechanisms at temperatures below — 37 °C 
(Liu and Penner, 2005). Although LP05 includes the im- 
mersion nucleation on soot it is not considered in this work 
as it is still highly uncertain (Karcher et al., 2007). In the 
LP05 parameterization, nucleated ice crystal number con- 
centration is computed based on temperature, relative hu- 
midity, sulfate and dust number concentration, and sub-grid 
updraft velocity. The sub-grid variability of updraft veloc- 
ity needed for driving the parameterization is derived from 


www.atmos-chem-phys.net/12/12061/2012/ 


Atmos. Chem. Phys., 12, 12061-12079, 2012 



12064 


the square root of turbulent kinetic energy (TKE) (Bretherton 
and Park, 2009) with an assumed maximum threshold value 
of 0.2 m s -1 (Gettelman et al., 2010). The limiting IN number 
concentration above which homogeneous nucleation is pre- 
vented is derived from parcel model simulations as a func- 
tion of temperature and sub-grid updraft velocity, and used 
in the parameterization of the competition between homoge- 
neous and heterogeneous ice nucleation. We find that the sub- 
grid variability of updraft velocity tends to be near 0.2 m s -1 
(the upper limit) much of the time in CAM5. Our analy- 
sis of vertical velocity measurements in the Interhemispheric 
Differences in Cirrus Properties From Anthropogenic Emis- 
sions (INCA) and Small Particles in Cirrus (SPartICus) cam- 
paigns indicates higher mean updraft velocities than the up- 
per threshold limit of 0.2 m s -1 (Zhang et al., 2012). This up- 
per threshold limit is also lower than that seen in the observa- 
tions in the tropical anvil cirrus of 0.30-0.5 m s -1 (Jensen et 
al., 2009). However, the model may predict too high sub-grid 
updraft velocity at high altitudes (e.g., near the tropopause) 
(Hoyle et al., 2005). Too high updraft velocity may produce 
anomalously high ice supersaturations, which may induce a 
homogeneous ice nucleation bias. The effects of sub-grid up- 
draft velocity on the ice nucleation in cirrus clouds will be 
investigated in our future study. 

In mixed-phase clouds with temperatures between 0 and 
—35 °C, contact freezing of cloud droplets through Brownian 
coagulation with insoluble IN (assumed to be mineral dust) 
(Young, 1974), deposition/condensation nucleation (Mey- 
ers et al., 1992), and immersion freezing of cloud droplets 
(Bigg, 1953) are included. Ice multiplication through Hallet- 
Mossop secondary ice production due to accretion of drops 
by snow is included following Cotton et al. (1986). More de- 
tailed description can be found in Gettelman et al. (2010). 

We note that there are a few changes between this study 
and Gettelman et al. (2010) in CAM5 for the LP05 parame- 
terization: (1) We use all Aitken mode sulfate aerosol for ho- 
mogeneous nucleation while Gettelman et al. (2010) used a 
portion of the Aitken mode sulfate with diameter greater than 
0.1 pm; (2) Accumulation mode dust (0.1 to 1.0 pm diame- 
ter) in addition to coarse mode dust (1.0-10.0 pm diameter) 
is used in the ice nucleation; and (3) We extend the LP05 pa- 
rameterization for ice nucleation in pure ice clouds with tem- 
peratures between 0 and — 37 °C where only heterogeneous 
nucleation occurs. This change is to be consistent with BN09 
parameterization which includes ice nucleation for pure ice 
clouds at temperatures less than — 37 °C as well as tempera- 
tures between 0 and — 37 °C. We note that the third modifi- 
cation has a small impact on model results (not shown). 

2.4 BN09 Parameterizations in CAM5 

For this study, we implemented the BN09 framework in 
CAM5. This framework explicitly considers the competition 
between homogeneous and heterogeneous freezing, and has 
the flexibility of using different IN “nucleation spectra” ei- 
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Table 1. CAM5 cases in this study. 


Case Name 

Description 

LP 

LP05, combined nucleation 

LPhom 

LP05, pure homogeneous nucleation 

LPhet 

LP05, pure heterogeneous nucleation 

BN 

BN09, combined nucleation 

BNhom 

BN09, pure homogeneous nucleation 

BNhet 

BN09, pure heterogeneous nucleation 


ther from empirical correlations derived from observation 
data or from CNT. The “nucleation spectrum” describes the 
number of ice crystals nucleated from IN as a function of ice 
supersaturation. In the version of the BN09 code used in this 
study, there are four options for “nucleation spectrum”. The 
first one (MY92) is derived from in-situ measurements of IN 
concentrations for temperature between 250 K and 266 K and 
ice supersaturation between 2 % and 25 % following Mey- 
ers et al. (1992). The second one (PDG07) is also based 
on in-situ measurement described by Phillips et al. (2007). 
The third one (PDA08), presented by Phillips et al. (2008), 
is a more comprehensive formulation that takes the surface 
area contribution from different types of aerosol (i.e., dust 
and soot) with sizes larger than 0.1 pm diameter and freez- 
ing modes (i.e., deposition and immersion) into account. The 
last one is derived from a semi-empirical form of CNT. The 
PDA08 is used in this study as it uses fewer tuning param- 
eters. Barahona et al. (2010) carried out sensitivity studies 
using a global chemical transport model with four different 
nucleation spectra and found that when the competition be- 
tween homogeneous and heterogeneous freezing is consid- 
ered, global mean ice crystal number concentration vary by 
up to a factor of twenty depending on the nucleation spec- 
trum used. 

Under the BN09 framework, the number of ice crystals 
nucleated is a function of aerosol number concentration, tem- 
perature, updraft velocity, and accommodation coefficient for 
water vapor deposition on ice crystals. Ice crystal number is 
derived from solving equations for temperature, ice super- 
saturation, and ice water content for an ascending air par- 
cel with initial condition that is saturated with respect to ice. 
For consistency, the same sub-grid updraft velocity is used to 
drive both BN09 and LP05. BN09 is applied only for pure ice 
clouds (i.e., cirrus clouds with temperatures less than — 37 °C 
as well as pure ice clouds with temperatures between 0 and 
— 37 °C when only heterogeneous nucleation takes place). In 
mixed-phase clouds, we use the same heterogeneous ice nu- 
cleation parameterizations as in Gettelman et al. (2010), as 
documented in Sect. 2.3. 
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3 Modeled ice crystal number concentrations 

We run CAM5 with prescribed present-day sea surface tem- 
peratures (SST) and aerosol and precursor emissions from 
Lamarque et al. (2010) for 5yr plus 3-month spin up, and 
results from the last 5-yr are used in the analysis. All sim- 
ulations have been carried out at 2.5° x 1.9° horizontal res- 
olution with 30 vertical levels. In the control case we use 
the LP05 and BN09 parameterizations in CAM5 including 
the competition between homogeneous and heterogeneous 
nucleation (i.e., “combined simulations”, denoted below as 
LP and BN, respectively). Two additional sets of simulations 
were performed: pure homogeneous freezing (hereafter as 
LPhom and BNhom, respectively), and pure heterogeneous 
freezing (hereafter as LPhet and BNhet, respectively), for 
LP05 and BN09 parameterizations by turning off their repre- 
sentation of heterogeneous or homogeneous nucleation, re- 
spectively. Table 1 gives the list of simulations carried out. 
We note that ice formation in mixed-phase clouds discussed 
in Sect. 2.3 keeps the same between all these simulations. 

The annual and zonal mean mass and number concentra- 
tions of mineral dust in the accumulation mode with the size 
(diameter) range of 0.1-1 pm and in the coarse mode with 
the size (diameter) range of 1-10 pm are shown in Fig. 1. 
About 90% of the total dust mass is in the coarse mode, 
while the accumulation mode dust dominates the total dust 
number concentration (Liu et al., 2012). Accumulation mode 
dust number concentration is in the range of 0. 1-1.0 cm -3 
in the NH upper troposphere (100-300 hPa), where coarse 
mode dust number concentration is one order of magnitude 
lower. Dust is uplifted from the NH subtropics and midlati- 
tude sources to the middle and upper troposphere and trans- 
ported towards the Arctic in the upper troposphere. The dust 
number concentration in the upper troposphere agrees rea- 
sonably well with the model result by Hendricks et al. (201 1). 
As we noted in Sect. 2, dust number concentration is used 
as an input for the heterogeneous nucleation formulations in 
LP05 and BN09. 

An offline (box model) comparison between the LP05 and 
BN09 parameterizations for the pure homogeneous nucle- 
ation case done by Barahona and Nenes (2008) showed that 
the two parameterizations give similar ice number concen- 
trations at cirrus temperatures less than — 40 °C and updraft 
velocity less than lms -1 . The comparison of in-cloud ice 
crystal number concentration (Ni) for LPhom and BNhom 
is shown in Fig. 2 (middle panel). We can see that distribu- 
tion patterns of Ni from LPhom and BNhom cases are sim- 
ilar, with higher Ni in the Southern Hemisphere (SH) and 
in the tropics. Both LPhom and BNhom cases predict Ni 
higher than 500 L -1 in the tropical tropopause and in SH. 
Ni is reduced at warmer temperatures in the middle tropo- 
sphere (300-500 hPa) with values of 20-100 L -1 , consistent 
with the nucleation rate reduction with increasing temper- 
atures for homogeneous ice nucleation (Koop et al., 2000; 
Pruppacher and Klett, 1997) and aggregational growth of 
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small ice crystals during the gravitational sedimentation from 
above. In the lower troposphere (700-850 hPa) Ni increases 
to 100-200 L -1 in storm track regions due to the ice nucle- 
ation in mixed-phase clouds. 

To compare BN09 and LP05 under the same condition, 
we turn off the soot contribution to IN number in BN09 (i.e., 
only dust particles can act as IN in both LP05 and BN09). 
An offline (box model) comparison showed that heteroge- 
neous IN number concentration from CNT-based spectrum is 
typically a few times higher than that from the PDA08 (Bara- 
hona and Nenes, 2009a) with the same input of aerosols. As 
shown in Fig. 2 (lower panel), Ni calculated from BNhet with 
PDA08 is in the range of 10-50 L -1 in the upper troposphere 
regime, much lower than that (100-300 L _1 ) in the NH mid- 
latitudes from LPhet. The higher Ni (100-300 L _1 ) in the 
subtropical and mid-latitudes of NH in LPhet is associated 
with the Saharan and Asian dust. The different Ni between 
LPhet and BNhet is due to the different formulations of IN 
number concentration (i.e., the CNT method in LPhet and 
the PD A0 8 in BNhet) even though they have the input of al- 
most the same dust number concentrations between the two 
simulations (not shown). 

Ice crystals formed from heterogeneous nucleation can in- 
hibit the subsequent homogeneous nucleation through water 
vapor consumption in the cloud parcel. Offline comparison 
by Barahona and Nenes (2009b) showed that the limiting het- 
erogeneous IN number (Mim) that would completely inhibit 
homogeneous nucleation is lower in LP05 parameterization 
than that in BN09 parameterization (i.e., homogeneous nu- 
cleation is more easily prevented in LP05 parameterization) 
under the same IN number concentration. Therefore we ex- 
pect that the occurrence frequency of homogeneous nucle- 
ation will be lower with the LP05 parameterization hav- 
ing both higher heterogeneous IN number concentration and 
lower Aii m than the BN09 parameterization when both ho- 
mogeneous and heterogeneous nucleation are allowed to oc- 
cur and compete. As a result Ni will be lower with LP05 
in regions where there is sufficient heterogeneous IN to in- 
hibit the homogeneous nucleation. Indeed, this is the case. 
As shown in Fig. 2 (upper panel), Ni in the NH upper tro- 
posphere is significantly reduced (by up to a factor of 2) in 
LP compared to LPhom due to higher heterogeneous dust IN 
concentration (100-300 L _1 ) there from LP05 (lower panel), 
clearly indicating the inhibition of homogeneous nucleation 
by heterogeneous nucleation. There are no obvious changes 
between BN and BNhom due to the lower heterogeneous dust 
IN number concentration (10-50 L -1 ) in the upper tropo- 
sphere which does not have a strong impact on the homo- 
geneous nucleation there, consistent with the study by Bara- 
hona et al. (2010). 

Figure 3 shows the frequency distribution of homoge- 
neous nucleation occurrence events in CAM5 (i.e., count as 
1 if homogeneous nucleation occurs, otherwise 0 at each 3- 
dimensional grid) for combined homogeneous and hetero- 
geneous nucleation case and pure homogeneous nucleation 
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Fig. 2. Comparison of in-cloud ice crystal number concentration (L -1 ) with LP05 (left) and BN09 (right) parameterizations for the combined 
homogeneous and heterogeneous nucleation (upper), pure homogeneous nucleation (middle), and pure heterogeneous nucleation (lower) 
simulations in CAM5. 
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Fig. 3. Frequency distribution of homogeneous nucleation occurrence events in CAM5 for combined homogeneous and heterogeneous 
nucleation case (upper) and pure homogeneous nucleation case (lower) with LP05 (left) and BN09 (right) parameterizations. 




Fig. 4. Percentage contribution from heterogeneous ice nucleation 
to total ice crystal number concentration diagnosed in the ice nucle- 
ation parameterizations for the combined LP and BN cases. 


case. We can see that the frequency of occurrence of homo- 
geneous nucleation is similar for LPhom and BNhom case 
(lower panel). Heterogeneous IN reduce the frequency of ho- 
mogeneous nucleation significantly for LP05 parameteriza- 
tion (upper panel), due to its higher heterogeneous IN num- 
ber concentration and lower 7Vii m , while the reduction is evi- 
dent only in the NH high-latitudes for BN09. 


Figure 4 shows the percentage contribution of ice crystals 
from heterogeneous ice nucleation to total ice number diag- 
nosed in the ice nucleation parameterizations for the com- 
bined LP and BN cases. We can see that homogeneous nu- 
cleation is the dominant contributor to Ni in the upper tro- 
posphere in the tropics and in the SH, and contribution from 
heterogeneous nucleation is much smaller there (< 10 %). In 
both LP and BN simulations, heterogeneous IN concentra- 
tion is not high enough to completely turn off the occur- 
rence of homogeneous freezing there. In the NH heteroge- 
neous nucleation becomes more important especially for the 
LP case. The contribution of heterogeneous nucleation to to- 
tal ice crystal number can reach over 50 % for the LP case 
due to the larger IN number from the CNT. The contribution 
from heterogeneous nucleation is less than 10 % (and mostly 
less than 1 %) in the upper troposphere for the BN case. Us- 
ing the probability distributions of RHi measured during the 
INCA campaign, Haag et al. (2003) suggested that in situ cir- 
rus clouds in the NH midlatitudes were likely formed from a 
combination of homogeneous and heterogeneous nucleation 
processes, whereas SH midlatitude cirrus clouds were dom- 
inated by the homogeneous nucleation. These INCA results 
appear to support the LP05 parameterization more than the 
BN09 parameterization with the heterogeneous IN spectrum 
of Phillips et al. (2008). We note that the relative importance 
of ice nucleation mechanisms depends on the in-cloud up- 
draft velocity. If there was a high updraft velocity bias in 
CAM5, reducing the bias would make it easier for IN to pre- 
vent RHi from reaching the threshold values for the initia- 
tion of homogeneous nucleation, and thus increase the im- 
portance of IN from PDA08 on the overall ice crystal num- 
ber. This will produce larger differences between BNhom 
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Fig. 5. In-cloud ice crystal number concentration (L -1 ) versus temperature. Model results are sampled every 3 hours over tropical, mid- 
latitude and Arctic regions including the observation locations reported in Kramer et al. (2009). The 50 % percentile (solid line), 25 % and 
75 % percentiles (error bar) are shown for each 1 K temperature bins. 


and BN results and improve the agreement of BN results with 
the INCA observations. 

4 Comparison of model results with observations 

Kramer et al. (2009) collected an extensive aircraft obser- 
vational dataset in the temperature range of 183-240 K. We 
compare model simulations against this dataset for pure ho- 
mogeneous, pure heterogeneous and combined cases with 
LP05 and BN09 parameterizations. We output modeled in- 
cloud Ni every 3 h and bin them into 1 K interval of temper- 
atures from 190 to 240 K (i.e., in the cirrus cloud regime) 
in the upper troposphere from 30° S to 75° N. The range be- 
tween 25th and 75th percentiles of modeled Ni within each 
temperature bin is shown in Fig. 5. The most distinctive fea- 
ture of Fig. 5 is that observed Ni tends to increase with in- 
creasing temperature for the whole temperature range. Ni is 
in the range of ~ 10-60L -1 at temperatures below 205 K. 
There have been studies to suggest that under these very 
low temperatures (i.e., near the tropical tropopause layer) 
aerosols rich with organic matter (e.g., Froyd et al., 2009) 
may become glassy, prevent homogeneous ice nucleation, 
and only allow heterogeneous nucleation to happen (Jensen 
et al., 2010; Murray et al., 2010). Barahona and Nenes (201 1) 
also suggested that low Ni would result from the emergence 
of dynamical equilibrium states, favored at low temperatures. 
There is no obvious Ni trend in the temperature range be- 
tween 205 K and 230 K, and Ni is higher than that in the 
lower temperatures and within the range of 20 L - 1 to a few 
hundred per liter. Above 230 K, Ni increases with tempera- 
ture and is frequently higher than 1000 L -1 . This could be 
due to the shattering of large ice crystals at warmer tempera- 
tures, as suggested in Field et al. (2003) and Lawson (201 1). 
Simulated Ni decreases with increasing temperature for pure 
homogeneous nucleation (LPhom and BNhom) and com- 
bined simulations (LP and BN) especially with BN09 pa- 


rameterization, consistent with the predominance of homo- 
geneous nucleation. Ni from LPhom is close to that from BN- 
hom. Ni from LP is lower than that from LPhom especially in 
the NH mid-latitudes with temperatures between 205-230 K 
by up to a factor of 2, because of the stronger effect of hetero- 
geneous ice nucleation in the LP05 parameterization than in 
BN09 (Fig. 4). Ni from BN is almost the same as BNhom due 
to the low heterogeneous IN concentration from the PDA08 
spectrum. Ni from pure homogeneous (LPhom and BNhom) 
and combined simulations (LP and BN) overestimates ob- 
served Ni by a factor of 10-100 in the temperature range 
of 190-205 K, while in the temperature range of 205-230 K, 
Ni from combined simulations is closer to the observed val- 
ues. Compared to pure homogeneous and combined nucle- 
ation simulations, pure heterogeneous nucleation simulation 
with the BN09 parameterization (BNhet) shows better agree- 
ment with the observation in the temperature range of 190- 
205 K, although it is not clear if the heterogeneous IN spec- 
trum (PDA08) used in BNhet can be extrapolated to such 
low temperature regime. Ni from pure heterogeneous nucle- 
ation simulations (LPhet and BNhet) is too low (20-40 L _1 ) 
in the temperature range of 205-230 K compared to the ob- 
servational data. LPhet which is based on CNT gives a higher 
Ni than BNhet which is based on PDA08 in the whole tem- 
perature range. The above comparison is consistent with the 
assertion that homogeneous ice nucleation is not able to ex- 
plain the very low Ni in the low temperature range of less 
than ~ 205 K near the tropical tropopause, and pure hetero- 
geneous nucleation may dominate there (Jensen et al., 2010). 
The low Ni in cirrus clouds at these temperatures may also be 
due to the occurrence of “dynamical equilibrium” states be- 
tween ice production and sedimentation loss (Barahona and 
Nenes, 2011). However, pure heterogeneous nucleation on 
dust with both LP05 and BN09 parameterization would not 
be able to explain observed Ni concentrations (~ 100 L _1 ) 
in the temperature range of 205-230 K, even in the LPhet 
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Fig. 6. Probability distribution frequency (PDF) of relative humidity in 150-300 hPa from 30° S to 60° N from pure homogeneous, pure 
heterogeneous and combined nucleation simulations for LP05 (left) and BN09 (right) parameterizations, in comparison with humidity ob- 
servations in the upper troposphere from the Measurements of Ozone, Water Vapor, Carbon Monoxide and Nitrogen Oxides by Airbus 
In-Service Aircraft (MOZAIC) and from the Atmospheric Infra-Red Sounder (AIRS) satellite. 


simulation with CNT. This indicates that homogeneous nu- 
cleation is often the dominant ice nucleation mechanism in 
this temperature range. Considering the large fluctuation of 
observed Ni at a given temperature in this temperature range, 
heterogeneous nucleation may still have important effects 
during many cloud events, especially at low cooling rates. 

Ice supersaturation is the driver of ice nucleation, and the 
level and distribution thereof gives insight on the relative im- 
portance of freezing mechanisms (e.g., Haag et al., 2003). 
We compare the probability distribution of relative humidity 
with respect to ice, RHi, in 150-300 hPa from 30° S to 60° N 
from instantaneous model output with pure homogeneous, 
pure heterogeneous and combined nucleation for LP05 and 
BN09 parameterizations, respectively. The modeled RHi is 
compared in Fig. 6 with humidity observations in the upper 
troposphere from the Measurements of Ozone, Water Vapor, 
Carbon Monoxide and Nitrogen Oxides by Airbus In-Service 
Aircraft (MOZAIC) project (Gierens et al., 1999) and from 
the Atmospheric Infra-Red Sounder (AIRS) satellite (Gettel- 
man et al., 2006a). Both observed probability distributions of 
RHi from MOZAIC and AIRS show exponential decreases 
with increasing of RH. Note that MOZAIC may be subject 
to compressional heating artifacts (inflating the RH), and 
may also be biased towards cloud-free conditions since pi- 
lots might avoid passages through thick cirrus layers. AIRS 
with a large RH uncertainty (~ 20 %) suffers from its 2-3 km 
vertical weighting function that may reduce the RH and the 
frequency of occurrence of thin supersaturated layers (Get- 
telman et al., 2006b). The all-sky model results can be bi- 
ased towards low supersaturations owing to the water vapor 
deposition on ice crystals to remove supersaturations inside 
clouds for CAM5 with a large model time step of 30 min. In- 
deed, the peak in frequency at 100 % RHi in the simulations 
is due to cloud points that are not sampled by MOZAIC or 
AIRS. Modeled supersaturation for pure homogeneous and 
combined nucleation with the BN09 parameterization (BN- 


hom and BN) is higher than AIRS and lower than MOZAIC. 
Results from BNhom and BN are similar. Model simulations 
with the LP05 parameterization (LPhom and LP) underes- 
timate the occurrence frequency of RHi larger than ~ 120- 
140%. As shown in Fig. 6, with the decrease of Ni from 
pure heterogeneous nucleation simulations (LPhet and BN- 
het), occurrence frequency of high ice supersaturation is sig- 
nificantly increased. The low Ni (20-50 L -1 ) in the upper 
troposphere from BNhet simulation could not consume water 
vapor efficiently by depositional growth of ice crystals and 
thus occurrence frequency of ice supers aturation from this 
simulation is larger than that from MOZIAC at RHi ~ 150- 
180 %. We note that low Ni (and thus small surface area of 
ice crystals) has been used to explain the persistent high su- 
persautration in cold cirrus clouds (Murray et al., 2010; Peter 
et al., 2006). 

Number concentration and size distribution of ice parti- 
cles were routinely measured over the US Department of 
Energy Atmospheric Radiation Measurement (ARM) South- 
ern Great Plains (SGP) site (36.6° N, 97.5° W) from De- 
cember 2009 to June 2010 during the SPartICus campaign 
(http://campaign.arm.gov/sparticus/). The suite of cloud par- 
ticle probes included a 2D-S probe (Lawson et al., 2006), 
which was equipped with probe tips that reduced shattering 
and a post-processing algorithm designed to remove residual 
shattered particles (Lawson, 2011). SPartICus provides rel- 
atively long-term statistics (~ 150 h of in situ data) of con- 
centrations of ice particles in cirrus. A recent comparison be- 
tween historical and SPartICus measurements in deep cirrus 
shows that SPartICus data typically contain lower concen- 
trations of small (< ~ 50 pm) ice particles, particularly near 
cloud base where large particles contribute the most to shat- 
tering (Lawson, 2011). Figure 7 shows the probability distri- 
bution of in-cloud ice number concentrations from SPartICus 
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Fig. 7. Probability distribution frequency (PDF) of in-cloud ice number concentration (L 1 ). Model results are sampled every 3h and 
interpolated into the field measurement site location (36.6° N, 97.5° W). All data are sampled over 7 km. 


data obtained during January to June 2010. 1 The most fre- 
quently observed Ni is in the range of about 10-200 L -1 , 
although ice number concentrations larger than 1000 L -1 are 
also observed in a few occasions. 

To compare with the aircraft measurements collected over 
the SGP area during the SPartICus campaign, we sampled in- 
stantaneous ice crystal number concentrations over the SGP 
site (interpolated using values in four GCM grid-boxes clos- 
est to the site) every three hours from January to June from 
5-yr simulations for different cases listed in Table 1. As most 
of the aircraft measurements were collected above 500 hPa, 
only samples above this level are considered in the analy- 
sis. All LP05 model simulations overestimate observed fre- 
quencies for Ni higher than ~ 100 L _1 , and underestimate 
observed frequencies for Ni less than 20-40 L _1 . This is also 
the case for BN and BNhom. BNhet produces two modes 
in the probability distribution: one at 10L -1 and another at 
200-400 L - 1 . The lower Ni mode is due to the heterogeneous 
ice nucleation, while the higher mode is due to the convec- 
tive detrainment, which also contributes to the occurrence of 
Ni at 200-400 L -1 for other cases (not shown). These com- 
parisons suggest that homogeneous nucleation may play an 
important role in cirrus clouds to produce the high frequency 
of Ni higher than 100 L -1 over the ARM SGP site during 
the SPartICus. However model simulations underestimate 
the frequencies of occurrence for Ni less than 40 L -1 ex- 
cept for the BNhet simulation. This may indicate that CAM5 
with large horizontal and vertical grid sizes may underes- 
timate the growth processes of ice crystals by aggregation 
and/or rimming, which can reduce the ice number concentra- 
tion in clouds (Spichtinger and Gierens, 2009a, b). Although 
LP and LPhet modeled histograms agree better with the ob- 
servations than that of LPhom, a conclusion on the dominant 
role of heterogeneous nucleation mode during the SPartICus 
is still premature due to the discrepancies between model 
simulations and observations. Future analysis will evaluate 

1 The SPartICus data used in this analysis are considered pre- 
liminary. However, final processing of the data is not expected to 
significantly alter the results shown here. 


the model representation of sub-grid updraft velocity (cool- 
ing rate), which is critical for the occurrence of homogeneous 
nucleation, and separate the in situ cirrus cases from the con- 
vective anvil cirrus cases in the comparison of model results 
and observations. 

5 Atmospheric impacts 

In this section, we will examine the impact of heterogeneous 
IN on upper tropospheric air temperature, specific humid- 
ity, and cloud properties by comparing simulations between 
the pure homogeneous nucleation and the combined homo- 
geneous and heterogeneous nucleation simulations. Figure 8 
shows latitude-pressure cross sections of annual mean Ni 
(grid-box average), IWC (grid-box average), cloud cover, 
temperature, and specific humidity from the LP simulation 
and differences of these variables between LPhom and LP. 
There is a statistically significant increase in Ni at the 90 % 
level of the Student’s t-test in LPhom compared to LP owing 
to the high number of heterogeneous IN in the NH midlati- 
tudes from dust sources in North African and Asian deserts 
(see Fig. 1). Compared to LP, IWC in the NH upper tro- 
posphere in LPhom is increased by 20% due to smaller 
ice crystal sizes and slower sedimentation from the ice- 
producing layers. The strong increase in Ni enhances cloud 
cover in the NH midlatitudes by 2-5 %. Tropical upper tro- 
pospheric temperature increases by ~ 0.5 K due to the long- 
wave radiative heating from the increased presence of cirrus. 
Upper tropospheric water vapor generally increases by 10- 
20%. The higher tropical tropopause temperature increases 
the water vapor saturation pressure, and this together with 
reduced sedimentation of cloud ice allows more water va- 
por to be transported to the lower stratosphere (Rosenfield 
et al., 1998). Thus specific humidity in the lower strato- 
sphere is increased by ~ 10%. Liu et al. (2009) and Hen- 
dricks et al. (2011) found similar impacts on the upper tro- 
pospheric and lower stratospheric (UTLS) water vapor mix- 
ing ratio when the ice crystal number concentration in cir- 
rus clouds is changed. In comparison, the difference in Ni 
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Fig. 8. Latitude-pressure cross sections of annual mean Ni (grid-box average), IWC (grid-box average), cloud cover, temperature, and 
specific humidity from combined ice nucleation simulation (left) and differences of these variables between pure homogeneous and combined 
nucleation simulations (right) for the LP05 parameterization in CAM5. Differences significant at the 90 % level of the Student’s t-test are 
depicted by dots. 


between BNhom and BN is much smaller with maximum 
changes in the tropical tropopause about 50 L -1 . Changes in 
IWC, cloud cover, temperature, and UTLS water vapor are 
statistically insignificant at the 90 % level between BN and 
BNhom (not shown). 


Figure 9 shows the annual and zonal mean shortwave 
cloud forcing (SWCF), longwave cloud forcing (LWCF), net 
cloud forcing (CF), high cloud cover (CLDHGH), ice wa- 
ter path (IWP), cloud- top ice effective radius (ACTREI), and 
column cloud ice number concentration (CDNUMI) from the 
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Table 2. Global annual mean cloud properties from sensitivity simulations listed in Table 1 for liquid water path (LWP, gm -2 ), ice water 
path (IWP, gm -2 ), shortwave cloud forcing (SWCF, W m -2 ), longwave cloud forcing (LWCF, W m“ 2 ), total cloud cover (CLDTOT, %), 
high cloud cover (CLDHGH, %), and column ice number concentration (CDNUMI, x 10 8 m“ 2 ). 


Case Name 

LWP 

IWP 

SWCF 

LWCF 

CLDTOT 

CLDHGH 

CDNUMI 

LP 

46.9 

18.3 

-53.0 

29.5 

65.4 

42.5 

4.82 

LPhom 

47.4 

19.4 

-55.0 

31.9 

65.6 

42.6 

5.86 

LPhet 

45.3 

17.2 

-46.3 

22.6 

62.1 

38.8 

1.34 

BN 

47.1 

17.8 

-52.9 

28.8 

64.5 

41.1 

5.75 

BNhom 

46.9 

17.8 

-52.9 

29.1 

64.7 

41.6 

5.79 

BNhet 

44.6 

15.4 

-43.8 

18.6 

61.0 

37.5 

0.67 


model simulations carried out in this study. Modeled SWCF, 
LWCF and CF are compared with satellite retrievals from 
the Clouds and the Earth’s Radiant Energy System (CERES) 
(Wielicki et al., 1996), and CLDHGH with combined Cloud- 
Sat and CALIPSO cloud fractions (Gettelman et al., 2010). 
Overall, SWCF from combined nucleation simulations (LP 
and BN) agrees well with the observations, except in the 
tropics and subtropics where the SWCF is too strong, proba- 
bly because of the too high ice crystal number concentration 
there in these simulations (Fig. 5). Noticeable differences in 
SWCF between different ice nucleation simulations are only 
in the tropics and subtropics. This is expected since SWCF 
in the extratropics is dominated from middle and low-level 
clouds. The SWCF differences reach ~ 10Wm -2 between 
different simulations for both LP05 and BN09 parameteriza- 
tions, with pure homogeneous nucleation simulations hav- 
ing the strongest (more negative) SWCF, followed by the 
slightly weaker SWCF from the combined nucleation simu- 
lations and then weakest SWCF from the pure heterogeneous 
nucleation simulations. This is consistent with the lowest Ni 
from pure heterogeneous nucleation simulations, the small- 
est IWP, and the weakest SWCF. Similar differences between 
different simulations occur for LWCF, but expand to the ex- 
tratropics, since cirrus clouds have an important contribu- 
tion to LWCF globally. LWCF from different model simu- 
lations underestimates the observations in the midlatitudes. 
The net CF from all the model simulations for both LP05 
and BN09 is stronger than observations in the tropics and 
subtropics. The net CF differences between different simula- 
tions are small due to the cancellation of SWCF and LWCF 
differences. The lowest Ni from the pure heterogeneous nu- 
cleation (LPhet and BNhet) results in the largest ice crystal 
size and highest sedimentation rate and thus the smallest IWP 
(Fig. 9b). The same is true for CLDHGH from BNhet which 
is lower by ~ 10 % than other two simulations (BN and BN- 
hom) in the SH extratropics due to the lowest Ni from BNhet 
simulation (Fig. 9a). There is no obvious difference in the 
cloud-top ice effective radius except for the BNhet simula- 
tion which has a larger ice effective radius due to the lower 
Ni (Fig. 9b). The cloud-top ice effective radius is similar 
among the three simulations for LP05, although differences 


are much larger at lower cloud levels. Column ice number 
concentration from the pure homogeneous nucleation simu- 
lation is higher than that from the combined nucleation sim- 
ulation especially in the NH midlatitudes for the LP05 pa- 
rameterization, and both simulations are a few times higher 
than that from the pure heterogeneous nucleation simulation 
for both LP05 and BN09 parameterizations. 

Table 2 gives the global and annual means of LWP, 
IWP, SWCF, LWCF, total and high cloud cover, and col- 
umn ice number concentrations from different simulations. 
Comparing combined (LP) and pure homogeneous nucle- 
ation (LPhom) simulations, heterogeneous dust IN reduces 
the column ice number concentration by 20 %. This reduces 
global mean SWCF (less negative) by 2.0 zb 0.1 W m -2 (1 a 
uncertainty) and LWCF (less positive) by 2.4 ±0.1W m -2 , 
respectively for LP05. The net change in global mean cloud 
forcing (CF) is -0.40 ± 0.20 Wm“ 2 (cooling). The SWCF 
and LWCF in the pure heterogeneous nucleation simulation 
(LPhet) are much weaker (by 6.7 ± 0.2 and 6.9 ± 0.2 W m -2 , 
respectively on the global mean) compared to the combined 
simulation, due to a factor of ~4 lower column ice num- 
ber concentration. The column ice number concentration 
with BNhet is half of that with LPhet. Comparing combined 
simulation (BN) and pure homogeneous nucleation (BN- 
hom), heterogeneous dust IN has a negligible effect on the 
global mean SWCF, and reduce the global mean LWCF by 
0.24±0.14Wm -2 (cooling). The SWCF and LWCF in the 
pure heterogeneous nucleation simulation (BNhet) are lower 
than the BN simulation by 9.1 ± 0.2 and 10.2 ±0.1W m -2 , 
respectively. 

Figure 10 shows the annual and zonal mean distributions 
of net CF differences between the combined and pure homo- 
geneous nucleation simulations (i.e., net CF changes due to 
the dust IN) and between the pure heterogeneous and com- 
bined nucleation simulations for LP05 and BN09. Table 3 
gives the net CF changes averaged over different latitude 
bands. Statistically significant changes caused by the dust IN 
occur in the NH midlatitudes (—1.81 ±0.47 Wm -2 forLP05 
and — 0.39 ±0.20 Wm -2 for BN09). The net CF changes 
between the pure heterogeneous and combined nucleation 
simulations are statistically significant in the tropics and in 
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Fig. 9a. Annual and zonal mean distributions of shortwave cloud forcing (SWCF), longwave cloud forcing (LWCF), net cloud forcing (CF), 
and high cloud cover (CLDHGH) from different model simulations with LP05 (left) and BN09 (right) parameterizations. Black-dashed lines 
are from observations. 


the SH for LP05 due to the low dust IN concentration there, 
while the net CF changes are significant at all latitude bands 
for BN09. 


6 Discussion and conclusions 

In this study the impact of heterogeneous dust IN on upper 
tropospheric cirrus clouds is investigated using CAM5 with 
two ice nucleation parameterization frameworks (LP05 and 
BN09). Both parameterizations consider the homogeneous 
and heterogeneous nucleation and competition between the 
two mechanisms. The heterogeneous IN concentration from 


LP05 is based on the CNT, and is higher than that from 
BN09, which is based on the PDA08 spectrum. The differ- 
ence is larger in the NH upper troposphere (by up to a factor 
of ten) where mineral dust particles from deserts are more 
abundant. Because of the inhibition of homogeneous nucle- 
ation by heterogeneous dust IN, the frequency occurrence of 
homogeneous nucleation in the upper troposphere is substan- 
tially reduced in the NH compared with that from the pure 
homogeneous nucleation simulation in LP05. The impact of 
heterogeneous dust IN is much smaller in BN09 because of 
the lower dust IN concentration as predicted by the PDA08. 
As a result, the overall ice crystal number concentration is re- 
duced by up to a factor of 2 in NH in the combined nucleation 
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Fig. 9b. Same as Fig. 9a but for ice water path (IWP), cloud-top ice effective radius (ACTREI), and column cloud ice number concentration 
(CDNUMI). 


Table 3. Global, tropical (30° S-30° N), Northern Hemisphere (NH) and Southern Hemisphere (SH), mid-latitude (30-60° S and 30-60° N), 
and high-latitude (60-90° S and 60-90° N) annual mean changes (A) of net cloud forcing (CF, W m -2 ) between the combined and pure 
homogeneous nucleation simulations and between the pure heterogeneous and combined nucleation simulations for the LP05 and BN09 
parameterizations. Standard deviations (±cr) are estimated from means of each of 5 yr. 


Case Name 

LP-LPhom 

LPhet-LP 

BN-BNhom 

BNhet-BN 

Global 

-0.40 ±0.20 

-0.17 ±0.28 

-0.24 ±0.28 

-1.22 ±0.08 

Tropics 

0.15 ±0.21 

0.78 ± 0.24 

-0.25 ±0.52 

0.30 ±0.19 

NH mid-latitudes 

-1.81 ±0.47 

0.10 ±0.29 

-0.39 ±0.20 

-2.29 ±0.28 

NH high-latitudes 

-1.41 ±1.35 

-0.54 ±0.60 

-0.14 ±0.47 

-2.69 ±0.59 

SH mid-latitudes 

-0.15 ±0.81 

-1.44 ±0.85 

-0.08 ±0.68 

-2.19 ±0.47 

SH high-latitudes 

-0.45 ±0.11 

-4.48 ±0.37 

-0.34 ±0.24 

-5.85 ±0.40 


simulation compared to the pure homogeneous nucleation 
simulation for LP05. Contribution of ice crystals nucleated 
from heterogeneous nucleation to total nucleated ice number 
concentration can reach over 50 % in the NH from the LP05 
parameterization compared to less than 10 % globally for the 
BN09 parameterization. 

We compare modeled cloud ice number concentration 
from different simulations with the two parameterizations 
with observations. Good agreement only occurs at tempera- 


tures lower than 205 K for the pure heterogeneous nucleation 
simulation with BN09 using the PDA08’s heterogeneous IN 
spectrum. However, pure heterogeneous nucleation simula- 
tions give too low ice number concentrations at the tem- 
perature range of 205-230 K, and thus homogeneous nucle- 
ation may play an important role in the ice formation. This 
is confirmed from the relatively long-term (~ 150 h) statis- 
tics of cirrus clouds obtained at a NH midlatitude site during 
the SPartICus campaign, which shows the large frequency 
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Fig. 10. Annual and zonal mean distributions of net cloud forcing (CF) differences between combined and pure homogeneous nucleation 
simulations (upper) and between pure heterogeneous and combined nucleation simulations (lower) with LP05 (left) and BN09 (right) param- 
eterizations. The vertical bars overlaid on solid lines depict the ranges of the two standard deviations (2a) estimated from the CF differences 
of each of 5 yr at different latitudes. 


(41 %) of high ice number concentrations (> 100 L" 1 ). How- 
ever, the large fluctuation (by more than one order of magni- 
tude) of observation data at a given temperature also suggest 
the role of heterogeneous nucleation affecting or completely 
inhibiting the homogeneous nucleation as well as the impor- 
tance of aggregation and rimming of ice crystals. Model sim- 
ulations with homogenous nucleation dominating the ice for- 
mation underestimate the occurrence frequency of ice num- 
ber concentration lower than 40 L -1 , while overestimate the 
frequency of ice number concentration higher than 100 L -1 
during the SPartICus. 

Heterogeneous dust IN may have a significant global 
impact on cloud radiative forcing (SWCF and LWCF) 
by 2Wm -2 , ice water content by 20%, and upper tro- 
pospheric water vapor mixing ratio by 10-20% if there 
is sufficient number of heterogeneous IN (~200L -1 ) in 
the upper troposphere. The net cloud forcing change is 
—0.40 zb 0.20 W m -2 . However, with a small number of het- 
erogeneous IN (20-50 L -1 ) in the upper troposphere the 
net cloud forcing change becomes statistically insignificant 
(— 0.24 ± 0.28 Wm -2 ). These results highlight the impor- 
tance of understanding the ice nucleation mechanisms and 
of quantifying the number concentrations and properties of 
heterogeneous IN. We note that the relative importance of 
heterogeneous IN on upper tropospheric cirrus clouds also 
depends on the cooling rate (i.e., vertical velocity) of air 
parcels in clouds. If the vertical velocity is sufficient high 
(e.g., larger than 0.5 ms -1 ) homogeneous nucleation may 


eventually take place after the heterogeneous nucleation and 
dominate the overall ice number concentration (e.g., Karcher 
and Lohmann, 2003). The evaluation of modeled sub-grid 
vertical velocity with SPartICus remote sensing and in situ 
aircraft observations will be conducted in our next study. 

Lohmann et al. (2008) also examined the competition ef- 
fects between heterogeneous nucleation from dust aerosol 
and homogeneous nucleation from sulfate aerosol using the 
ECHAM5 model. They assumed that the homogeneous nu- 
cleation will be completely switched off (i.e., only heteroge- 
neous nucleation can occur) when IN concentration is higher 
than a threshold value of 1L -1 , and homogeneous nucle- 
ation occurs elsewhere. By this they found a much stronger 
net CF cooling effect (— 2.0 Wm -2 ) due to dust IN than 
— 0.40 ± 0.20 Wm -2 obtained in our study. In our study, 
when the IN concentration is about 20 L -1 (as calculated 
from PDA08 in BNhet shown in Fig. 2), homogeneous nu- 
cleation still plays a dominant role in the combined simu- 
lation for the BN09 parameterization (Fig. 4). Only when 
the IN concentration is higher than 100 L -1 in the North- 
ern Hemisphere (as calculated in LPhet shown in Fig. 2), the 
heterogeneous nucleation then becomes very important in the 
combined simulation for the LP05 parameterization (Fig. 4). 
The much lower IN threshold value assumed in Lohmann 
et al. (2008) to completely switch off the homogeneous nu- 
cleation can be the reason for the larger net CF cooling ef- 
fects of dust IN found in Lohmann et al. (2008). This is con- 
firmed from the significant reduction of ice crystal number 
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concentration on the global scale in their combined simula- 
tion as compared to their pure homogeneous nucleation sim- 
ulation and the similar ice crystal number in their combined 
simulation to that in their pure heterogeneous nucleation sim- 
ulation as shown in Fig. 2 of Lohmann et al. (2008). The 
global annual mean vertically integrated ice crystal number 
is reduced by 47 % in their combined simulation as compared 
to their pure homogeneous nucleation simulation (Table 2 of 
Lohmann et al., 2008), while in our study the reduction is 
only 1 8 % for the LP05 parameterization and the reduction 
is only seen in the Northern Hemisphere (Fig. 2 and Table 2 
of this study). This indicates a much greater (and a domi- 
nant) role of heterogeneous nucleation in the combined sim- 
ulation of Lohmann et al. (2008). A comparison of the sub- 
grid cloud updraft velocity between Lohmann et al. (2008) 
and this study will help to understand the sensitivity of ice 
nucleation to aerosols, which will be done in our future study. 

Model evaluation with observations indicates several ar- 
eas of future model improvement. The sub-grid features of 
ice supersaturation, and ice nucleation and growth need to be 
better represented in global climate models with horizontal 
grid sizes of 100 km. These ice microphysics parameteriza- 
tions need to be linked with the cloud macrophysics scheme 
in climate models. This can be achieved by implementing 
a statistical cloud macrophysics scheme (e.g., Karcher and 
Burkhardt, 2008) into CAM3 (Wang and Penner, 2010) and 
CAM5 (M. Wang, personal communication, 2012), which 
will be evaluated with observation data. Further model im- 
provement would also separate the prediction of cloud ice 
crystals by sources, such as the ice crystals formed from ho- 
mogeneous nucleation, heterogeneous nucleation and from 
convective detrainment (Hendricks et al., 2011). BN09 ice 
nucleation parameterization will also be modified to account 
for the impact of pre-existing ice crystals on ice nucleation 
and to include ice nucleation on glassy aerosols (Murray et 
al., 2010) and on ammonium sulfate (Abbatt et al., 2006). As 
we find from this study, major differences in CAM5 simula- 
tions with LP05 (the default scheme in CAM5) and the new 
implemented BN09 are due to their different heterogeneous 
ice nucleation formulations. We note that BN09 has four IN 
spectra to choose. The comparison of these IN spectra on ice 
crystal number concentration and impact on cirrus clouds un- 
der the same BN09 framework will be presented in our future 
study. 
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